Potential Energy Surfaces of Actinide Nuclei from a Multi-dimensional Constraint 
Covariant Density Functional Theory: Barrier Heights and Saddle Point Shapes 
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For the first time the potential energy surfaces of actinide nuclei in the (/32o, /322, /^ao) deformation 
space are obtained from a multi-dimensional constrained covariant density functional theory. With 
this newly developed theory we are able to explore the importance of the triaxial and octupole shapes 
simultaneously along the whole fission path. It is found that besides the octupole deformation, the 
triaxiality also plays an important role upon the second fission barriers. The outer barrier as well as 
the inner barrier are lowered by the triaxial deformation compared with axially symmetric results. 
This lowering effect for the reflection asymmetric outer barrier is 0.5 ~ 1 MeV, accounting for 
10 ~ 20% of the barrier height. With the inclusion of the triaxial deformation, a good agreement 
with the data for the outer barriers of actinide nuclei is achieved. 

PACS numbers: 21.60.Jz, 24.75. +i, 25.85.-w, 27.90.-fb 



Since the first interpretation of the nuclear fission by 
the barrier penetration [Q, it has been a difficult task 
to describe this phenomenon theoretically. For example, 
in order to study the fission problem, one should ffi^st 
have very accurate information about the fission barrier; 
a 1 MeV difference in the fission barrier could result in 
several orders of magnitude difference in the fission half- 
life. Particularly, to explore the island of stability of 
superheavy nuclei (SHN) |||-^, it is more and more de- 
sirable to have accurate predictions of fission barriers of 
SHN ff. 

Nowadays three types of models are used for calculat- 
ing fission barriers. During a long period the majority 
of these works is based on the macroscopic-microscopic 
(MM) models |, | |ll|. The MM models make use of the 
Strutinsky shell correction method, allowing fast calcula- 
tions of multi-dimensional PES's containing most of the 
important shape degrees of freedom. Until now it is still 
an important candidate for large scale fission barrier cal- 
culations based on the examination of multi-dimensional 
PES's 1 10 . In recent years, the rapid development of the 



density functional theories (DFT) also makes it possible 
to calculate the fission barriers fully self-consistently 
[Tst . There are mainly two reasons to apply DFT's in the 
study of fission properties. First, many new functional 
forms and effective interactions are proposed with much 
better performances for the excited state as well as the 
ground state calculations [p^|-pO[ . Fission barrier calcula- 
tions are also helpful for developing these DFT's. Second, 
in DFT's much more shape degrees of freedom can be 
included self-consistently. For example, the symmetry- 
unrestricted Skyrme-Hartree-Fock-Bogoliubov model has 
been applied for the fission studies |2^. Besides these 
two types of models, there exist also methods intending 
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to combine the advantages of the MM and self-consistent 
models, such as the extended Thomas- Fermi method [ p2[ . 

The double-humped fission barriers of actinide nuclei 
can be used to benchmark the predictive power of the- 
oretical models ||lj, p3|-p5t. Various shape degrees of 
freedom play important and different roles in the occur- 
rence and in determining the heights of the inner and 
outer barriers. For examples, it has long been known 
from MM model calculations that the inner fission bar- 
rier is usually lowered when the triaxial deformation is 
allowed, while for the outer barrier the reflection asym- 
metric (RA) shape is favored |2^. Later on these 
points were also revealed in the non-relativistic and 
relativistic ^ density functional calculations, re- 
spectively. It is thus customary to consider only the tri- 
axial and reflection symmetric (RS) shapes for the in- 
ner barrier and axially symmetric and RA shapes for 
the outer one [|[ |2^, It has been pointed out that 
"there is no reason for a fissioning actinide nucleus not 
to penetrate all symmetry-breaking shapes on its way 
from the first (triaxial) to the second (mass-asymmetric) 
saddle" The non-axial octupole deformations are 

considered in both the MM models and the non- 
relativistic Hartree-Fock theories However, a multi- 
dimensional structure of PES's including both the triax- 
ial and RA shape degrees of freedom has not been ex- 
plored yet in the framework of covariant DFT. In this 
paper, we will investigate the influence of the triaxiality 
and the octupole shape on the PES's all the way from 
the ground state to the flssion configuration when both 
shape degrees of freedom are included simultaneously. To 
this end, not only as many as self-consistent symmetries 
should be broken, but also multi-dimensional constraints 
arc needed [ p4[ . 

To calculate the potential energy surfaces and fission 
barriers, in this work we use the covariant density func- 
tional theory (CDFT) |l7[ |l|, H, |§ . By breaking 
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not only the axial Q but also the reflection sym- 
metries |39t] , we developed multi-dimensional constraint 
CDFT's in which the functional can be one of the follow- 
ing four forms: the meson exchange or point-coupling 
nucleon interactions combined with the non-linear or 
density-dependent couplings [Q. If not specified, the 
functional form of the point-coupling nucleon interaction 
with non-linear self energy terms and the parameter set 
PC-PKl are used in this work. 

For the parametrization of the nuclear shape, we adopt 
the conventional ansatz in mean field calculations. 
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where Qa^i is the mass multipole operators. When the 
axial and reflection symmetries are broken, the nuclear 
shape is invariant under the reversion of x and y axes. 
In other words, the intrinsic symmetry group is V4 and 
all shape degrees of freedom /3x^ with fi even, includ- 
ing the triaxial (/j, 7^ 0) and octupole (A = 3) deforma- 
tions, are possible. Irrespective with the self-consistent 
symmetries, the single particle wave functions and vari- 
ous densities are expanded on an axially deformed har- 
monic oscillator basis I 
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In order to get a fast 
convergence of the results against the basis size, in the 
elongated direction more states are included in the ba- 
sis. Following Warda et al. |4^, the basis is truncated as: 
nz/Qz + i'2,n± + \m\)/Qj_ < Ncut- Here , nj_ , and m are 
quantum numbers characterizing each state in the basis, 
Qz = MAX(l,6,/6o), Qx = MAX(1, &i/&o), and 60, 6„ 
and b± are the oscillator lengthes. The calculated bind- 
ing energy of ^^°Pu at /32o = 1-3 varies only about 130 
keV and 20 keV when A'cut increases from 16 to 18 and 
from 18 to 20. This means a good convergence and such 
a truncation scheme with A'cut = 16 ensures a 0.2 MeV 
accuracy for the deformation range we are interested in. 
In the present work, A^'cut = 16 (20) is used in the triax- 
ial (axial) calculations. More details of the convergence 
study will be given in Rcf. The BCS approach is im- 
plemented in our model to take into account the pairing 
effect. Since it has been found that the BCS calculation 
with a constant pairing gap can not provide an adequate 
description of the fission barriers , we use a delta force 
for the pairing interaction with a smooth cutoff 41, H. 



We performed one- (1-d), two- (2-d), and three- 
dimensional (3-d) constrained calculations for the ac- 
tinide nucleus ^"^^Pu. In Fig. |l| we show the 1-d poten- 
tial energy curves (PEC) from an oblate shape with /32o 
about —0.2 to the fission configuration with /32o beyond 
2.0 which are obtained from calculations with different 
self-consistent symmetries imposed: the axial (AS) or 
triaxial (TS) symmetries combined with reffcction sym- 
metric (RS) or asymmetric cases. The importance of 
the triaxial deformation on the inner barrier and that of 
the octupole deformation on the outer barrier stressed 

kS] arc clearly seen here: The 
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triaxial deformation reduces the inner barrier height by 



FIG. 1. (Color online) Potential energy curves of 2*°Pu with 
various self-consistent symmetries imposed. The solid black 
curve represents the calculated fission path with V4 symmetry 
imposed, the red dashed curve that with axial symmetry (AS) 
imposed, the green dotted curve that with reflection symme- 
try (RS) imposed, the violet dot-dashed line that with both 
symmetries (AS & RS) imposed. The empirical inner (outer) 
barrier height Bcmp is denoted by the grey square (circle). 
The energy is normalized with respect to the binding energy 
of the ground state. The parameter set used is PC-PKl. 



more than 2 MeV and results in a better agreement with 
the empirical datum; the RA shape is favored beyond 
the fission isomer and lowers very much the outer fis- 
sion barrier. Besides these features, we observe for the 
first time that the outer barrier is also considerably low- 
ered by about 1 MeV when the triaxial deformation is 
allowed. Again, a better reproduction of the empirical 
barrier height can be seen for the outer barrier. We note 
that this feature can only be found when the axial and 
reflection symmetries are simultaneously broken. 

How the PES of ^*°Pu becomes unstable against the 
triaxial distortion can be seen much more clear in Fig. |^ 
in which we show 2-d PES's from calculations without 
and with the triaxial deformation. When the triaxial 
deformation is allowed, the binding energy of ^*°Pu as- 
sumes its lowest possible value at each (/320j/33o) point. 
At some points we get non-zero /322 values. That is, non- 
axial solutions are favored at these points than the axial 
ones. The triaxial deformation appears mainly in two 
regions in Fig. ^j. One region starts from the first saddle 
point and extends roughly along the direction of the (3^0 
axis up to a very asymmetric shape with /^ao ^ 1.0. In 
this region the values of /322 are about 0.06 ~ 0.12 cor- 
responding to 7 ^ 10°. The energy, especially the inner 
barrier height, is lowered by about 2 MeV. The other re- 
gion is around the outer barrier and the /322 values are 
about 0.02 ~ 0.03 corresponding to 7 ^--^2°. About 1 
MeV is gained for the binding energy at the second sad- 
dle point due to the triaxiality. In other regions, e.g., in 
the ground state and fission isomer valleys, only axially 
symmetric solutions arc obtained. 

Next we examine the full 3-d PES of ^'*"Pu obtained 
from the newly developed multi-dimensional constraint 
CDFT. For simplicity, in Fig. | arc shown only five typ- 
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FIG. 2. (Color online) Potential energy surfaces of Pu in 
the (;02O, (Sao) plane from calculations (a) without and (b) with 
the triaxial deformation included. The energy is normalized 
with respect to the binding energy of the ground state. The 
numbers in (b) show the values of /322 at these points. The 
fission path is represented by a dash-dotted line. The ground 
state and fission isomer are denoted by full and open circles. 
The first and second saddle points are denoted by full and 
open triangles. The contour interval is 1 MeV. 
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FIG. 3. (Color online) Sections of the three-dimensional PES 
p£240p^ in the (/322, /Sao) plane calculated at /32o= 0.3 (around 
the ground state), 0.6 (around the first saddle point), 0.9 
(around the fission isomer), 1.3 (around the second saddle 
point) and 1.6 (beyond the outer barrier), respectively. The 
energy is normalized with respect to the binding energy of the 
ground state. The contour interval is 0.5 MeV. Local minima 
are denoted by crosses. 



ical sections of the 3-d PES of 240pu in the (/322,;S3o) 
plane calculated at /320= 0.3 (around the ground state), 
0.6 (around the first saddle point), 0.9 (around the fission 
isomer), 1.3 (around the second saddle point) and 1.6 (be- 
yond the outer barrier), respectively. Many conclusions 
can be drawn by examining these 3-d PES's. First, the 
ground state and the fission isomer arc both axially and 
reflection symmetric as what is shown in the 1-d PEC and 
the 2-d PES. But with the 3-d PES one can investigate 
the stabiUty of ^^''Pu against the P22 and /Sso deforma- 
tions. One finds that the stiffness of the fission isomer is 
much larger than that of the ground state against both 
the /322 and /33o deformations. Second, while around the 
inner barrier the shape of ^""^"Pu is triaxial and reflec- 
tion symmetric, the second saddle point which is close to 
P20 = 1.3 appears as both triaxial and reflection asym- 
metric shape. Third, the triaxial distortion appears only 
on the top of the fission barriers. 

It has been pointed out that one may obtain spurious 
saddle points if only a small number of shape degrees 
of freedom are constrained, see, e.g., Ref. Q. That is, 
the calculated fission path may jump from one valley 
to another and results in discontinuities in the lower- 
dimensional PES's; in some cases, a continuous path 
may even cross a higher saddle point. Although the 
spurious saddle points may not be excluded completely, 
most of them can be avoided if (1) the obtained fission 
path keeps to be continuous in the energy as well as the 
most important shape degrees of freedom and (2) the 
results are examined by higher-dimensional calculations. 
We have carefully checked the full 3-d PES and found 
that the fission path enters and exits the triaxial con- 
figuration rather smoothly, which tells that no sudden 
jump is found and the 1-d (with the /32o deformation 
constrained and /322,/330 deformations imposed) and 2- 
d (with /320j/?30 deformations constrained and the /322 
deformation imposed) calculations of the fission barriers 
may be well justified for ^^°Pu. It is clear that the con- 
tinuity of the fission path found in a lower-dimensional 
constraint calculation is a necessary but not sufficient 
condition for locating the correct saddle point. In order 
to have a strictly definite conclusion, one certainly should 
carry out multi-dimensional constraint calculations with 
even highcr-multipolarity deformations included. 

For the RS calculations, the triaxiality also lowers the 
fission path by a few MeV beyond the second saddle 
point. This point is illustrated by the dotted line in Fig. |l| 
and the local minima with /Sso = 0.0 in the /32o = 1-6 sub- 
figures of Fig. ^. However, it is relatively unimportant, 
because the RA fission is still the most favored one even 
when triaxiality is included. 

Guided by the features found in the 1-d, 2-d, and 3-d 
PES's of ^^°Pu, the fission barrier heights are extracted 
for even-even actinide nuclei whose empirical values are 
recommended in RIPL-3 (see Table XI in |Q ) . The em- 
phasis is put on the influence of the triaxial deformation 
on the two fission barriers. 

As it has been shown previously, around the inner bar- 
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FIG. 4. (Color online) The inner (Bf) and outer (Bf) barrier 
heights of even-even actinide nuclei. The axial (triaxial) re- 
sults are denoted by open (full) symbols. The empirical values 
are taken from Ref. pa] and represented by grey squares. 



rier an actinide nucleus assumes triaxial and reflection 
symmetric shapes. Thus in order to obtain the inner fis- 
sion barrier height we can safely make a one-dimensional 
constraint calculation with the triaxial deformation al- 
lowed and the reflection symmetry imposed. In Fig. ||(a) 
we present the calculated inner barrier heights and 
compare them with the empirical values. It is seen that 
the triaxiality lowers the inner barrier heights of these 
actinide nuclei by 1 4 MeV as what has been shown 
in Rcf . . In general the agreement of our calculation 
results with the empirical ones is very good with excep- 
tions in the two thorium isotopes and ^"^^U. For ^^°Th and 
^^^Th, the calculated inner barrier heights are smaller by 
about 2 or 1 MeV than the empirical values depending 
on whether the triaxial deformation is allowed or not. In 
these two nuclei, the outer barrier is higher than the inner 
one. This may result in some uncertainties when deter- 
mining empirically the height of the inner barrier which 
is not the primary one [ [l2[ . Similar results for ^"^"Th 
and ^•^^Th were obtained from the Skyrme-Hartree-Fock- 
Bogoliubov model [ p^ and very small inner barrier height 
was got for ^^^Th in Ref. For ^ssu, Sj from the ax- 
ial calculation agrees with the empirical value very well. 
The triaxiality reduces the barrier height by about 1.5 
MeV, thus bringing a discrepancy which was similar to 
the result in Ref. KB. 

To obtain the outer fission barrier height Bf, the sit- 
uation becomes more complicated because more shape 
degrees of freedom have important influences around the 
outer fission barrier. For example, the inclusion of the 
refiection asymmetric shape makes it possible to have in 
the (/320, Pso) plane two or more competing fission paths 
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FIG. 5. (Color online) (a) One-dimensional potential energy 
curve E ~ /320 and (b) two-dimensional potential energy sur- 
face E ~ (/?20, /33o) of ^''^Cm the outer barrier region with the 
axial symmetry imposed in the calculation. In both figures, 
the energy is normalized with respect to the binding energy 
of the ground state. The fission path I/II is represented by 
full/dash dot lines and the corresponding saddle point is de- 
noted by up/down triangles. The fission isomer is denoted 
by open circle. In (a), the empirical outer barrier height is 
depicted by dash-dotted line. In (b), the contour interval is 
0.5 MeV. 



with different octupole deformations. In consequence, 
one often observes in the 1-d PEC two or more fission 

paths. This happens in 244,246p^ 244,246,248^^^ 

the present study and we present a typical example in 
Fig. I for 2'**Cm. In this figure, one finds that there are 
two fission paths both in the 1-d E ~ /32o curve and in 
the 2-d E - (/32o,,33o) PES. One path denoted by "I" 
favors shapes with larger octupole deformations and the 
other denoted by "11" favors less RA shapes. In such 
nuclei, it is not safe to perform a 1-d constraint calcula- 
tion in order to get Bf. Thus we first assume the axial 
symmetry and make a 2-d calculation in the (/32o,/33o) 
plane from which we can approximately identify the low- 
est fission path /?3o™'''(/52o) and the location of the sec- 
ond saddle point. Then along this fission path, we per- 
form a 1-d /320-constraint calculation with the triaxial 
and octupole deformations allowed. At each point with 
1320, the initial deformations are taken as /^j"' = and 
^ini- /3^7'="'(/32o)- In this 1-d PEC, we can locate the 
second saddle point and extract the outer barrier height 
for each nucleus. 

In the lower panel of Fig. ^ we show the results of outer 
barrier heights Bf and compare them with empirical val- 
ues. For most of the nuclei investigated here, the triaxi- 
ality lowers the outer barrier by 0.5 ^ 1 MeV, accounting 
for about 10 ~ 20% of the barrier height. One finds that 
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our calculation with the triaxiality agrees well with the 
empirical values and the only exception is ^^®Cm. From 
the calculation with the axial symmetry imposed, the 
outer barrier height of ^*®Cm is already smaller than the 
empirical value. The reason for this discrepancy may be 
related to that there are two possible fission paths beyond 
the first barrier, as seen in Fig. ^. For the path I with 
a lower saddle point from which we get the outer fission 
barrier height, the barrier is very wide and for the path II 
with a higher saddle point, the barrier is relatively nar- 
row. Therefore the empirical value of the outer fission 
barrier height may not be easily extracted for the follow- 
ing two reasons: (i) There must be a strong competition 
between the two fission paths; (ii) When the empirical 
value of the outer barrier height is evaluated, it is usually 
assumed that the second barrier is in an anti-parabolic 
shape with a fixed and smaller width p5| . 

We also examined the parameter dependency of our 
results. The lowering effect of the triaxiality on the outer 
fission barrier is also observed when parameter sets other 
than PC-PKl are used. 

In summary, a multi-dimensional constrained covariant 



density functional theory is developed which allows us to 
study the importance of the triaxial and octupole shapes 
simultaneously along the whole fission path. The one- 
dimensional PEC E ~ /320, two-dimensional PES E ^ 
(;^20, ko), and three-dimensional PES E ~ l322,l33o) 
of actinide nuclei are shown and studied in details. Both 
the triaxiality and the reflection asymmetry plays cru- 
cial roles at and around the second saddle point. The 
outer barrier as well as the inner barrier are lowered by 
the triaxial deformation compared with axially symmet- 
ric results. For most of the nuclei investigated here, the 
triaxiality lowers the outer barrier by 0.5 ~ 1 MeV, ac- 
counting for about 10 ^ 20% of the barrier height. The 
calculated results of the outer barrier heights agree well 
with the empirical values. 
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